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Magnetic Domain Interactions 

By R. L. GRAHAM 

(Manuscript received March 18, 1970) 

In this paper, we initiate a study into the combinatorial aspects of a 
model of the interactions between discrete magnetic domains and their 
potential use in information processing devices. Starting with a simple 
model suggested by W. Shockley, ice demonstrate certain {surprising) 
capabilities as well as inherent limitations upon the possible applications 
of the interactions described by this model. It should be noted that this simple 
model does not take into account all of the possible interactions between 
magnetic domains. 

I. INTRODUCTION 

The subject of discrete magnetic domains in certain orthoferrite 
materials has been under active investigation during the past several 
years, both from a theoretical physical viewpoint as well as that of the 
device-oriented physicist (for example, see Refs. 1-6). Considerable 
progress has resulted from these efforts, although needless to say, the 
end is certainly not in sight. Particular attention has been directed 
toward the problem of applying this new technology to the very im- 
portant area of information processing devices, an area in which it 
seems to have natural and significant applications. 1,7 It is our intention 
in this paper to examine certain mathematical aspects of these applica- 
tions for a simple model of magnetic domain interactions suggested by 
W. Shockley. 

II. DESCRIPTION OF THE MODEL 

We shall begin by giving a very brief description of the physical 
situation and its translation into the mathematical model under con- 
sideration. The reader whose interests motivate him to seek a more 
technical explanation is referred to Refs. 6 or 8. 

Roughly speaking, thin platelets of certain orthoferrite materials 
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possess the property that under suitable (magnetic) conditions, small 
(~ 3 mils) discrete cylindrical magnetic domains, hereafter called 
"bubbles", may be stably supported. Moreover, these bubbles may be 
manipulated by the application of external magnetic fields as well as 
by their own mutual interaction (which in general causes two bubbles 
to repel one another). In a suitable physical environment, the location 
of a bubble in a piece of orthoferrite can be restricted to a finite set of 
possible positions within the material; these are ordinarily arranged in 
a rectangular array. It is possible to apply a local magnetic field to 
specific locations within the array with the following results: 1 

(i) If a bubble already occupies the position at which the field was 

applied, then nothing happens. 
(it) If no bubble occupies the position at which the field was applied 
and no bubble occupies any "nearby" position as well, then 
(still) nothing happens. 
(iii) If no bubble occupies the position at which the field was applied 
but at least one bubble occupies some "nearby" position, then 
some bubble at a nearby position will leave its original position 
and now occupy the position selected by the field. 

To eliminate the annoying indeterminancy in item (Hi) it is possible 
to apply "holding" fields to all but one of the "nearby" sites which has 
the effect that only a bubble at the unheld position can move. 

The mathematical model which will correspond to the preceding 
description will be phrased in the terminology of graph theory. The 
discrete positions at which bubbles may lie correspond to the set V of 
vertices of a graph G. Two sites which are "nearby" or "adjacent" to 
one another (this is assumed to be a symmetric relation) correspond to 
two vertices of G which are joined by an edge of G. Suppose bubbles are 
located at (the sites corresponding to) the subset X of vertices V. We 
define a command to be a directed edge e = (v 1 , v 2 ) with Vi and v 2 adjacent 
vertices of G. The command e transforms the locations of the bubbles 
from X to X e where 

v : \X - M U {*} if Ml, v,4X; 



{X, otherwise. 
In other words, if there is a bubble at v y but no bubble at v 2 and the 



t Of course, "careless" application of a magnetic field to an orthoferrite with 
bubbles can annihilate bubbles, create bubbles, split bubbles in two, deform bubbles 
into strips, and so on; but these pathological (though certainly useful) operations 
will not be considered in our model. 
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command e = (Vi , v 2 ) is applied to X, then the bubble at v t is moved 
to v 2 . Otherwise, the command e has no effect on X. A program is 
defined to be sequence P = (e t , e 2 , • • • , e,) of commands e, . In 
general, a program P maps the set 2 of all subsets of V into itself 
by X p = (• • -(X' 1 )") • • •)"■ It is the purpose of this paper to inves- 
tigate the mathematical properties of these maps. 

III. SOME BASIC PROPERTIES OF PROGRAMS 

We begin by making the assumption that G is the complete graph on 
n vertices, that is, all pairs of vertices of G are joined by an edge. 1. As 
mentioned in the previous section, a program P is a sequence of directed 
edges (e, , e 2 , • • • , e r ) and P acts on a subset X of the vertices V of G by 

x p = (••• ((x e y)---r 

where for e = (u, v'), 

(Z - {v} U (»') if vtX, v' iX; 



X' = 

[X, otherwise. 

If X C V then \X\ denotes the cardinality of X. We note 
Fact 1: For all X C 7, and all programs P, |X P | = |X| . 
This follows immediately from the definition of X p . 
The first interesting result we state is due to W. Shockley who 

called it the 

Non-decreasing Overlap Theorem: (Shockley) For all X, Y C V and all 

programs P, 

I x P n y* i ^ i x n y i ■ 

Proo/: Assume for some P = (c, , ■ • ■ , e r ) and subsets X, F C 7 we 
have|X P n 5' P | < I * D F | . Since X p = (•• -((Z*)") ■••)", there 
must exist a Zeasf j such that 

| x n+* p| yi-,.. i < i X i-i p, Y Pi | 

where P t denotes the program (?i , ••• , e A ). Thus, for X = X p ', 
Y = Y Pi and e = e, + r = (a, b) we have 

1 1' H Y' | < | X fl f I • 



t Nothing essential is lost by this simplifying assumption. The vertices and edges 
of the present model should not be confused with any incidental physical vertices 
or edges in a particular device. An edge of the model may be generated for example 
by transferring bubbles from a storage zone to an interaction zone and then returning 
the resultant to the storage zone. 
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Ifc?*a,c^b then ctX f) $ implies c e 1° f) Y\ If either a t £ f) ? 
ovbzl f\ t but not both then & e £' f] ?'• If both atJtfl ? and 
6 e j£ f) ? tnen « e ^° D ^ e ana " b*X e f) f\ Hence, in any case 

I f n f I * I * n 1 1 

which is a contradiction. This proves the theorem. 

Shockley noted that this result shows that there is no replicating 
program P*. By a replicating program, we mean the following: Starting 
with two fixed sets of vertices V and V" with V C\ V" = 
and 1-to-l map 0: V" — > F', we require that for each X C F, 

X p< fl V = X H V and 0(X P * fl V") = X C\V. 

In other words, P* does not disturb X f) F' and in F", P* creates a 
"copy" of X F'. 

To show this, suppose there were such a program P*. By choosing 
two subsets X and X' differing in a single element of F, their images 
X p ' and X' 1 " must differ in two points, namely, one in V and the 
corresponding point (under 0) in V". This, however, contradicts the 
non-decreasing overlap (NDO) theorem and therefore P* cannot exist. 

Another consequence of the NDO theorem is the nonexistence of a 
program P + which performs binary addition in the following way. 

Suppose V denotes a set of m ^ 1 pairs of vertices of G, V" denotes 
another set of m pairs of vertices disjoint from V, and V" denotes a 
set of m + 1 pairs of vertices, disjoint from V and V". We can imagine 
these sets arranged as shown in Fig. 1. 

We can represent an integer M, ^ M < 2 m , in the m pairs of V 
by letting the ;/'th pair of V denote the ;th binary digit in the binary 
expansion of M. This can be done, for example, by assuming that 

O U 
for each pair , either U £ X, U x i X, which will correspond to 

U U i 

a 0, or U Q $ X, Ui t X, which will correspond to a 1 . Thus, for m = 5 
the configuration [V x , U s , U z , F 4 , F 5 } (Fig. 2) would denote the 
integer 1001 1 (2) - 19. 

The addition program P + would operate by starting with V" in some 
fixed configuration (for example, all zeros) and with arbitrary integers 
M', M" loaded into V, V", respectively, to form the initial state X; 
after applying P* to X we should get the sum M' + M " in V". 

The reason that P + cannot exist as described is precisely that the 
NDO theorem would be violated. For consider the two additions: 
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Fig. 1 — Symbolic arrangement of vertex locations for addition. 

+ (2'" — 1) = 2'" — 1 and 1 + (2"' — 1) = 2'". The initial configurations 
differ in only two positions. The final configurations differ in at least 



m + 1 however, since 2'" — 1 = 11 • • • 1 (2) and 2'" = 100 • ■ • (2) . 
Thus, by the NDO theorem we get a contradiction and our assertion 
is proved. 

We give another example of a program which does not exist. If 
e = (a, b) is a command and a, b e X then X e = X. In the case that 
a and b are both in X, we say that there is interference as e acts on A'. 
(We can think of the bubble at b as interfering with the attempted 
movement of the bubble at a to vertex b.) Similarly, if P = (ei , • • • , e n ) 
we say that there is interference as P acts on X if for some i there is 
interference as e { acts on X" "'. We note 

Fad 2: If P acts on X with no interference then 

X p = U \x] P . 

ztX 

Proof: It is sufficient to establish this for the case P = e = (a, b). 
In this case 
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Fig. 2 — A typical configuration representing an integer. 
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16, if x = a; 



)X) = 

[x, otherwise. 
Thus 



X - {a} U {6}, ifaeX; 



U [x} p = 

" x [X, otherwise. 

But by the hypothesis of no interference, we cannot have both a and 
6 t X. Thus 

x p = \X~ [a] U {&}, if atX = y {x]P 
[X, otherwise 

and the fact is established. 

FactS: For X = {a, b, c, z } , there does not exist a program P such that 

{a, b} p = ( 0| «}, 

{b,c} p = {a,z\, 

{c,a} p = {&,*}. 

Prooj: Suppose such a P exists. If P acts on these sets with no 
interference then we would have by Fact 2, 

{c,z} = {{a} p ,{b} P \, 

{a,z\ = \{b) P Ac\ P \, 

{Ml = {\c} P ,ia} P \, 

which is impossible since the union of the left-hand sides of the equations 
cannot equal the union of the right-hand sides. Thus, if P = (d , • • • , e„) 
we may assume that there is a least i, 1 ^ i ^ n, with P<_! = (ei , 
e<_i) such that e t acts on at least one of the sets {a, b} p, ~ l , {6, c\ 
{c, a} Pi ~ l with interference. To be specific, assume that it is the set 
{a, &}*'"*, that is, e> = ({a} Pi ~ x , {b} Pi ~ l ) (the other two cases are 
similar) . By Fact 2 we have 

(0,6}'"'- {{0}"-, {&}"-}. 

(«,a}"-« He}"", {a)"-}. 

Therefore 

{b,c} Pt = ({&}**-.{«}**-}"- (W"-\M""'I 



ifi-l 
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and 

\c,a\ Pi = [{e}"-\ (o}""}" = i\c\ Pi -\ (&}"-}■ 
Hence, 

\a,z) = {b,c\ P = [c,a} p = \b,z\ 

which is a contradiction. This proves the Fact 3. 

Note that the nonexistence of the program of Fact 3 does not follow 
directly from Fact 1 or the NDO theorem. A similar argument can be 
given to show that for X = (a, b, c,d,A, B, C, D, z) there is no program 
P such that 

\a,c} P = {A,z\, 
{a,d\ p = \B,z), 
\b,c\ p = {C,z}, 
\b,d\ p = {D,z\. 

IV. THE 2-VALUED BOOLEAN FUNCTIONS 

Our attention will now be focussed on the positive aspects of the 
model. In particular we shall be concerned with the problem of repre- 
senting the Boolean functions of m variables with appropriate programs. 
The way in which a function is to be represented is as follows. Suppose 
m = 2 and consider the function /: {0, 1} X {0, 1} — > {0, 1} by 

x y j(x, y) 
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If the values 1 and are interpreted as "true" and "false", respectively, 
then / is just the truth function of the familiar operation of alternation. 
V will be the set of six vertices (x ,x l ,y ,yi,fo, /i) which we indicate 
in Fig. 3. It is not difficult to show that no generality is lost by assuming 
there are no additional vertices. In fact, by using the pair of positions 
x , Xi in which to observe the result of the program, instead of providing 
the separate positions / , /i , it is true that if a Boolean function of 
m ^ 2 variables can be represented by a program in this general way, 
then it can be represented using just 2m vertices. The program P(f) 
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Fig. 3 — Symbolic arrangement of vertex locations for computing Boolean functions 

of two variables. 

which represents / is required to have the property that 

/„e{x, ,y } Pi/ \ Ui\x x ,y Q ] PU \ 

The correspondence between the indices of the vertices of V and the 
values of the variables of / is immediate. In terms of bubbles, one may 
think of the configurations shown in Fig. 4 as representing a and 1 
respectively (compare Fig. 2); P(j) is required to map each of the four 
possible initial states of the x,-pair and y,-pair into the correct value 
in the /,-pair. 

It is not difficult in this case to find an appropriate P(f), for example, 
we can take 

P (/) = Oo , ?/o)Oo , fo)(x l , ?/i)(v/i , /i). 

This is easily checked, as shown in Table I. We can write the preceding 
result in the shorthand form 



/ 



P(1) 



(0, 0, 0, 1) (x , y )(xo , /o)(xi , 2/i)(z/i , U) • 

Note that if /is defined by J'(x, y) = 1 — j{x, y), that is, /is the com- 
plement of /, then we can take 

P(f) = P(f)(x , X 1 )(X , tf ,)(«o , ?/ 2 )(/i , Xo)(fo , fi)(Xo , /o) 



o 



--BUBBLE--, 



o 



Fig. 4 — Configurations which represent and 1. 
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Table I— C 


umulative Effect 


OF P(/) 






(s, y) 


(xo, J/o) 


(•To, /o) 


(*ii ?/i) 


(»l,/l) 




/(*, v) 


(0,0) *-> \x ,y \ 

(0,1) ~ \x ,yi\ 

(1.0) «-» [x u y a \ 

(1.1) « |x„ J/,| 


(So, J/ol 
I.Vi, .'/"I 
|x , J/ol 
(Xli ?/il 


l/o, 2/o| 

h/ii yd 

fxi, J/ol 

l«lj 2/1 1 


l/o, yd 

1 2/i, ,'/o) 
l?/i, J/ol 
(si, 2/i 1 


l/o, ?/o| 

l/l, I/O 

1/., 2/o 1 


< ► 

< • 




i 
i 
i 



as a program which represents / (we leave this to the reader to verify). 
Table II, together with this remark about /, show that all of the 
16 possible 2-valued Boolean functions of two variables can be repre- 
sented by programs. 

A question which naturally arises at this point is whether all Boolean 
functions of m variables can be represented by programs in this manner. 
For m = 1, the answer is in the affirmative (the specific programs are 
left to the reader to discover); for m = 2, we have given the required 
16 programs; for m = 3, the answer is in the affirmative but the number 
(2 21 = 256) of programs prohibits their listing here; for m = 4, the 
answer is once again in the affirmative but the calculations necessary to 
establish this are much too long to be exhibited (there are, after all, 
2 24 = 65536 functions to consider). The cases m = 3 and m = 4 were 
established by J. H. Spencer." 

One may note that since all Boolean functions of two variables can 
be represented, then in particular the Sheffer stroke function given by 

x y Kr, y) 
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Table II — Programs for Boolean Functions of 2 Variables 



/ 


P(f) 


K), 0, 0,0) 


(xu, /o) (a: i, /o) 


(0, 0, 0, 1) 


(xi, 2/i ) (xi, /i) (x„, /o) (2/0, /o) 


(0, o, i, o) 


(xi, yo) (xi, /i) (x , /o) (2/i, /o) 


(0, 0, 1, 1) 


(x,, /i) (Xo, /n) 


(0, 1, 0, 0) 


(lo, ?/i ) (x , /i) (xi, /o) (2/o, /o) 


(0, 1,0, 1) 


(2/0, /o) (2/1, /i) ,. 


(0, 1, 1,0) 


(xo, yo) (xo, /o) (xi, ?/i) (2/1, 2/o) (xi, /o) (2/i, /i) 


(0, 1, 1, 1) 


(x , J/o) (x , /o) (xi, /i) (2/1, /i) 
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can also be represented. It is well known that any Boolean function 
of m variables can be generated by expressions containing just the 
variables and the stroke function. 10 Hence, one is tempted to conclude 
that any Boolean function is representable by a program. The flaw in 
this line of reasoning is that in order to express a particular Boolean 
function in terms of the stroke function, many occurrences of the stroke 
function and the variables are usually required. This in turn requires 
many "copies" of the variables to be available to the program in order 
to represent /. But we initially have only one pair of positions which 
indicates the value of any particular variable and by the NDO theorem 
we have seen that there cannot exist a "replication" program which 
would form extra copies of the values of the variables. Hence, within 
this model, we cannot use this technique to generate all the Boolean 
functions. It is certainly true however that if the model were extended 
to include bubble interactions which would allow replication of con- 
figurations (and such are known to exist physically), then all Boolean 
functions of ra variables could be represented exactly in the manner 
described. 

These initial results create considerable optimism concerning the 
possibility of representing all the Boolean functions of m variables. 
Such hopes are shattered however by the result (which we later prove) 
that there exists a Boolean function of 11 variables -which cannot be re-pre- 
sented by any program of this type. In fact, even though the fraction of 
the total number of Boolean functions of 11 variables which can be 
represented by programs can be shown to be < 10 -163 , the author is 
currently unable to exhibit any specific function which cannot be 
represented. Clearly, our understanding of this is less than complete. 
It is not unreasonable to hope that the representable functions could 
eventually be effectively characterized. 

We now restrict ourselves (without loss of generality) to representing 
the Boolean functions of m variables in the following way. We shall 
take V = {.r-j , x[ , x 2 , x' 2 , • • • x m , x„) to be a set of 2m vertices which 
we imagine to be arranged in pairs as illustrated in Fig. 5. As before, a 
bubble in the £,(£,') location of the pair (#,- , .r') will denote that the ith 
variable of the function / has the value 0(1). The way in which a program 
P(J) represents / is as follows. Choose a distinguished vertex atV. There 
is an obvious 1-1 correspondence between {0, l} m and the class C of 
all subsets X C V such that X intersects each {x t , x'i) in exactly one 
element given by 

a — (fli , • • • , a,„) <-> \y ( t V : y { = x { if a, = 0, 

y t - x{ if a, - 1, 1 si 1 £ m) = X. 
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Fig. 5 — Symbolic arrangement of vertex locations for computing Boolean functions 

of m variables. 



Let A< C {0, 1}", i = 0, 1, be the set of all a z {0, 1 } m such that /(a) = i, 
and let C, be the corresponding subsets of C. Our object is to find a 
program P(j) which distinguishes between the sets C Q and C, . (Note 
that Co U C t = C.) Specifically we shall say that P(j) represents / if 

at A 



for all 
for all 



XtC , 
XtC, . 



a $ X 

Let C denote the subset of all subsets x Q V with | x \ = in and 
for x and y distinct elements of V, let C(x) be the set of elements of C 
which contain x with C{y) denned similarly. 1 Consider the effect of 
the command {x, y) on the members of C(x) and C(y). There are four 
cases: 



C(y). 



(i) X c C(x), X e C(y). 

Then X ( " 0> = X and X ( " -) * C(x), X l *' u) 
(ii) XtC{x),XiC(y). 

ThenZ' 1 "" - X - [x\ U M and X to ' rt ^(i),!" 1 " e Cfo). 
(m) Z * C(x), X e C%). 

Then X lx ' v) = X and Z' 1 "" t C(x), X l '" ] e C(y). 
(w) XtC(x),X4C(y). 

Then Z*-- 1 - Z and I (I '"' * C(z), X* M 4 C(y). 

Hence, after the application of (x, y) to all the sets in C, the new sets 
C'(x), C'(y) (which now consist of all the subsets in C which contain x 
and y respectively) are related to C(x) and C(y) by 

C'(x) = C(x) fl C(y), 

C'fo) = C(x) U C(2/). 

Stated in these terms, the object of the program P(/) is finally to have 
C'"' (a) dC = C'o after it has been applied to all the sets in C. 

We give an example which illustrates these concepts. Let / be the 
Boolean function of three variables defined by: 



This approach was first suggested by J. H. Spencer. 
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x y z f(x, y, z) 
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V = [xi , x[ , x a , xi , x 3 , xi] and we take a = x[ . 

Co = { [Xt , X a , X 3 ] , [Xi , x' 2 , X a }, [Xj , x'a , Xa] , [x[ , X, , X z ] } , 

Cy = [{Xt, x 2 ,x' a ], {xi ,x 2 ,x^\, {xi t xi t ak} t [xi t xi, x' a ) ) . 
A program P(J) which achieves the separation is 

P(f) = (xl , x 2 )(x{ , X 3 )(X2 , Xz)(Xi , Xz)(X! , xi). 
That is, 

XeCo=*« = x[ tX Pin , 

A i Co => X[ f A 

If C(x) denotes the initial subset of C consisting of all the sets in C which 
contain x then we may conveniently record the sequential changes 
which occur in each current C(x) in terms of the original C(y)'s as 
the successive commands of P(f) are applied as shown in Table III. 
A little computation shows that the final set in the xj-row, the final 
C(x[), when intersected with C gives exactly 

{ [Xi , X 2 , Xa], [Xi , x' 2 , X a ], {Xi , X 2 , Xs\ , {xi ,X 2 , X 3 ] ] 

which equals C as required. 

In general the problem of representing Boolean functions reduces to 
the following problem. We start with the 2m classes C (0) (y) — C(y) (~) C, 
y t V. We are then allowed to replace two of the classes C (0) (y) and 
C m (y') by two (possibly) new classes C (0) (y) f| C w (y') and C m {y) \J 
C w (y'). We can repeat this operation as many times as desired with 
any pair of classes currently in the list. Our objective is to eventually 
generate a specified subset C* of C. 

We have already mentioned that for in = 1, 2, 3 and 4 it is possible 
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Table III — Cumulative Effect of P(f) 
(x,', x.) Or/, x z ) 



X\\ 


C(xi) 


x,': 


CW) 


x 2 : 


C(x s ) 


x 2 ': 


cW) 


x 3 : 


C(x 3 ) 


x 3 ': 


cW) 



C(x.) C(x.) 

C(x,') n C(x 2 ) C(x,') n C(x 2 ) n C(x 3 ) 

CW) u C(x 2 ) CW) u C(x«) 

CW) CW) 

C(x 3 ) (CW) n C(x 2 )) u CW) 

CW) CW) 

W, X 3 ) (Xi, x 3 ) 



x,: C'(xi) ('(lOn^i.Ou^di'JnC^jJJuCd 

x/: C(x l ')r.C(x 2 )nC'(x 3 ) C(x,')nC(x 2 ) nC(x 3 ) 

x ? : CW)uCW C(xx')uC(x 2 ) 

C(x 2 ')n((C(x,')r.C(x 2 ))uC(x 3 )) C(x 2 ')n ((CW) nC(i,)) uC(x 3 )) 

C'(x 2 ') u (CW) n C'(x 2 )) u C(x 3 ) CW u C*(x 2 ') u (CW) n C(x 2 )) u C(: 



*. 



- ...... . (x 8 ) 

C(x 3 ') CW) 

(x,, X,') 



x,: C(x.) n (CW) u (CW) n C(x 2 )) u C(x 3 )) n C(x,') n C(x 2 ) r> C(x 3 ) 

x,': (C(xi) n (C(x 2 ') u (CW) n C(x 2 ))) u C(x 3 )) u (C(x,') n C(x 2 ) n C(x 3 )) 

x 2 : C(xi') u C(x 2 ) 

x 2 ': CW) n ((CW) n C(x 2 )) u C(x 3 )) 

x 3 : C(xO u C(x 2 ') u (C(x/) n C(x 2 )) u C(x 3 ) 

x 3 ': C(x,') 



to generate a/*?/ subset of C in this manner. We proceed to show that 
for m = 11, there is a subset of C which cannot be generated. We first 
need several preliminary observations. 

To begin with, for a,b i V, let A and B denote the current sets C (,) (a) 
and C U) (b), respectively, after the ith command of the program P has 
been executed. In. other words, at this point in time C (a) is the class 
of all the original subsets of C which now contain a. For example, if 
a = x! z in the preceding example, then after the fifth (and final) com- 
mand of P(f), C (5) (x 2 ) is C w (x' 2 ) (] (C m (x' 2 ) f| C w (x 2 ) \J C M (x 3 )). It is 
immediate that if C U) (a) C C l,) (6) then the application of the command 
(a, b) as the (i + l)-st command of the program changes nothing. Hence 
we can assume that we only use commands (a, b) for ivhich at the time of 
their application C li) (a) £ C u \b) $ C i0 (a) (we say that C U) (a) and 
C (,, (6) are incomparable). 

Initially all the starting classes C m (x), x e V, are mutually incom- 
parable. In general suppose we have a family of classes D = {A,- ; 
1 ^ i ' ^ t), Ai C C, with exactly r of the ( 2 ) pairs of A, being comparable 
and assume A l and A 2 are incomparable. Consider the family D' = 
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D - {Ax} - [A 2 ] U {Ax f]A 1 } U {Ax U ^2}. We wish to determine 
how many pairs of the classes of D' are comparable. By definition 
D' = {Ax f) A 2 , Ai (J A 2 , A 3 , A A , • -, A t ). Of course for i, j ^ 3, 
the comparability between the classes A { and A t remains unchanged. 
There are several cases: 



(i) A< 3ii, ,A t 3 A 2 . 

Then A, 3 A, \J A 2 , A,- ^ A t f) A 2 
(it) 4,- 2 Ax , A, 3. A 2 . 

Then A, 3 A t fl ^2 
(m) Ai Ji.^iDi,. 

Then A, 3 Ai fl ^2 
(it;) A,- C Ax , Ai C A 2 . 

Then A, C Ax fl ^2 , 4* £ A t (J A. 
(v) Ai C Ax , Ai £ A 2 . 

Then A, Ci,|JA 2 
(vi) Ai $ Ax , A, C A 2 . 

Then A, C A t (J ^2 

Finally, we have a most important new comparability in D', namely 
Ax O A 2 Q Ax \J A 2 . Thus, at least r + 1 pairs of classes of D' are 
comparable. An immediate consequence of this observation is 

Fact 4-' We can assume that no program P(f) consists of more than ( 2 2 ) 
commands. 

Proof: Since after i (nontrivial) commands of a program P(f) have 
been applied, we must have (by induction) at least i pairs of the classes 
C <x) (x), x c V, being comparable and since there are just 2m classes and 
therefore ( 2 2 ) pairs of classes, then P(j) must have ^ ( 2 2 ) commands. 

Theorem. There exists a Boolean function of 11 variables which cannot 
be represented by a program. 

Proof: It is sufficient to show that for m = 11, there is a subset C* 
of C which cannot be generated by starting with the 2m classes C w (x), 
x e V, and recursively applying the transformation A, B — » A (~) B, 
A U B. Consider a typical program P = (ex , e 2 , • • • , e t ) and the corre- 
sponding expressions C U) (t), presented in Table IV. 

In choosing the ith command e,- of P there are at most ( z ™) — i + 1 
possibilities for e,- since after (ei , • • • , e.-O has been applied, at least i — 1 
of the pairs C''" 11 ^), C l *~ v (y) are comparable and thus neither (x, y) 
nor (y, x) can be the next command ei . Therefore there are at most 



CB r. 

n 



[( 2 2 m ) - *• + l] = [m(2m-l)]\ 
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Table IV — Cumulative Effect of P 

P: e\ e 2 e, e t 

«i: CW(xi) C«>(z,) C< 2 >(xi) ••• C<*>(si) C<*Hxi) 

xi'\ C™W) C« 1 >(xi') C«>(zi') ••• C»6d') C»>(zi') 

z m : C<°'(x m ) C<»(z m ) C™(x m ) ■■■ dHxm) C«KXm) 

x m ': C(°'W) C<»(x„') C< 2 >(x„,') ••• C«W) C"'(x m ') 



choices for the sequence of e, , since t ^ ( 2 ™) = m(2m — 1) by Fact 4. 
Also, for i ^ 1, each column C <0 (x), ieF, contains at most two new 
classes which did not occur in the preceding column since only two 
classes are changed at each step. Hence there are at most 

[m(2m- 1)] l( 2 ^) + 2m 

classes which can be generated by these rules where the additional 
term +2m comes from the 2m initial sets C {0) (x), x c V. On the other 
hand, since C contains 2 m sets X C V, then there are 2 2 " subsets of C 
which we must try to generate. We are doomed to failure however since 

l[m(2m - 1)]! (^™) + 2m| / ' 2 2 " -> 

as m — > oo . We list these expressions for several small values of m in 
Table V. Thus, not only are we guaranteed a single Boolean function 
of 11 variables which cannot be represented by a program, but in fact 
we have at least 10 015 of them. It seems quite likely that there exist 
Boolean functions of five variables which cannot be represented. 
However, at present, no specific example of a Boolean function is known 
which cannot be represented by a program. 



Table V — Bounds on the Number of Boolean Functions which 

Can Be Generated 



m 



[m(2m -l)]-'( 2 ^) +2?n 2= m 



2 
3 

10 
11 



4324 16 

19615115520006 256 

>10 3S8 <10 309 

<10 453 >10 616 
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V. SOME REMARKS 

A number of partial results are known concerning the preceding 
problems which we shall only mention briefly here. 

The generation of Boolean functions as described has the following 
very natural geometrical interpretation. For a fixed integer n, consider 
the set of the 2" vertices of an n dimensional cube C" and let A x , ■ ■ • , A 2n 
represent the 2n sets of 2 n_1 vertices which each lie on one (n — 1)- 
dimensional "face". In other words, if the vertices of C" are labelled by 
binary n-tuples in the usual way, then each A { corresponds to a set of 
2" -1 n-tuples in which some component is constant. As before, we are 
allowed to replace any tw^o sets A and B in the class of 2n sets of A P) B 
and A \J B. We can repeat this transformation as often as desired. The 
question is: which subsets ICC" can be generated in this manner. We 
have shown that there exists a set X C C 11 which cannot be so generated. 

More generally, suppose we start with a class of n formal sets 
X x , • • • , X n and ask which formal expressions in the X, can be generated 
using the transformation X, Y — > X (~| Y, X \J Y iteratively. It can be 
shown 11 for example, that all the elementary symmetric functions (using 
n and U in place of the usual • and +) can be generated. Let us call 
a w r ell-formed expression E in the X k 's symmetric in X< and Xj if the 
substitution X, — > X, , X, — > X, , leaves E unchanged. Thus we can 
write E in the form 

e = (x,- n x< n w,) u ((x f u Xi) n wj u w 3 

where the W { are well-formed (possibly empty) expressions in the X k 's 
not involving X,- or X,- . We say that we collapse X { and Xj in E if we 
apply the transformation X,- f] X,- — > X, , X,- (J X,- — > X,- , to form 

w = (i,. n w x ) u (x, n w 2 ) u w 3 . 

Certainly, if E can be generated using the transformations X, 
Y — > X n Y)X\J Y starting from Xj , • • • , X n , then there is a sequence 
of collapses starting with E and ending with some single variable X< . 
A basic theorem can be proved which asserts that if it is possible to 
generate E, then no matter how w r e collapse symmetric variables 
starting with the expression E we must reach some single variable X,- . 
In other words in attempting to collapse E to a single variable, we can 
never make a "bad" move. Once the structure of the expressions E 
which can be generated is sufficiently well understood, perhaps the 
representable subsets of C n can then be determined. 

Another line of research suggested by this bubble model is in the 
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following direction. For binary sequences x and y, define d(x, y), the 
(Hamming) distance between x and y, to be the number of positions in 
which the sequences x and y differ. The fact which prevented the 
existence of a program which could add two integers expressed to the 
base 2 was the fact that there are pairs of additions in which the binary 
expansions of the corresponding summands are close together (in the 
metric d) but whose sums are not close, thus conflicting with the NDO 
theorem. What we would like is a mapping m — * r(m) of integers into 
binary sequences for which we have 

d(r(m), r(n)) + d( T (m'), r(n')) ^ dQr(m + n), r (m' + n')). 

With only this constraint there are trivial solutions, for example, 

TO-+111 ••• 1. 

With this mapping we are essentially expressing m to the base 1 (well- 
known by many cultures to be inefficient for representing large numbers, 
say, those exceeding 10). Hence, we might require in addition that the 
number of binary sequences of length t which are in the range of the 
mapping r to be at least a for some fixed a > 1. Is it possible to find 
a suitable t for which an addition program is possible in this model of 
bubble interactions? 

Finally, we have just considered just one rather simple model in this 
paper. Physically, many other bubble interactions are possible (although 
some presently operate with significantly smaller margins than others) 
and this of course would lead to other models. It would be very interest- 
ing to understand the corresponding questions in some of these other 
models. 
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